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COLLISION ORBITS FOR A HILL'S TYPE PROBLEM 



CESAR CASTILHO 



Abstract. We study the planar problem of two satellites attracted by a center 
C " ;) ' of force. Assuming that the center of mass of the two-satellite system is on a 

C " ;) ' circular orbit around the center of force and using Levi-Civita regularization 

^^ ^ we prove the existence of an almost periodic orbit with an infinite number of 

^_^ collision between the satellites. 

^O ' Key Words: Hill's Problem, Regularization, Collisions. 

1. Introduction 

^ ■ Consider the planar problem of three bodies of masses mo,TOi and m2 in the case 

'nJ" I where mi and 1712 are much smaller than mg. The mutual atraction of the two small 

bodies can be usually neglected and the problem reduces in a fair approximation 

to two independent two-body problems. However if the distance between the two 

(^ • small bodies is small their mutual atraction can no longer be ignorated. This is 

known as Hill's problem [H . The derivation of Hill's equations usually found in the 
^D ■ literature assumes a hierarchy of masses for the three bodies: 



Too >> mi >> 7712, 



I ' and proceeds in two steps: First the limit m2 — > is taken, which gives the restricted 

'■^ I three-body problem; then the limit toi — > is taken. 

C^ ' A number of problems in celestial mechanics can be approximated by Hill's 

equations. Examples are: the Sun-Earth-Moon problem and the interaction be- 
tween satellites on nearby orbits. The nearby orbits problem was our motivation 
for this work. From the data obtained during the passage of Voyager near Saturn, 
j^ ' it was found that the two satellites Epimatheus and Janus around Saturn have al- 

H I most circular orbits, with radii 151,422 Km and 151,472 Km, and periods 16.664 

and 16.672. At each close encounter, they interchange orbits. To explain this be- 
haviour some models have been proposed, among them Scaling Techniques |CH| 



and Assymptotic Expansions [HP|. A common feature of the proposed models is 
an assymptotic expansion in the variable that represents the distance between the 
two satellites and truncation of higher order terms. In other words, such models 
assume the hyphothesis that there is a smaller bound for the distance among the 
satellites. In this paper we show that this hypotesis can not be assumed without 
further assumptions. Using Levi-Civita regularization we show the existence of an 
almost periodic orbit with an infinite number of collisions among the satellites. 
Therefore, it's not true in general that such lower bound exist. 

We take a different approach to Hill's problem. First, instead of taking the limits 
TO-i, m2 -^ we fix the body of mass tuq at the origin and assume mg = 1. We take 
TO-i, 7712 << 1 but no hierarchy for the masses mi and to,2 will be assumed. Second, 
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we assume that the center of mass of the two-satelhte system is on a circular orbit 
around the center of force. This is equivalent to consider the circular Hill problem 
ffl . The collision between the satellites is them regularized using the canonical form 



of the Levi-Civita regularization |SS]. The organization of the paper is as follows: 
In section 2 we develop our model for the Hill's Problem. The model is regu- 
larized and contrary to the usual study of Hill's problem we do not use a rotating 
system of coordinates (synodical). We show that the potential is symmetric with 
respect to reflections about the origin in the regularized physical plane. This will 
be crucial to prove existence of the collision orbit. 

In section 3 we study the collinear equilibria predicted by our model. It's an 
interesting feature of the model that the regularized equations allows to find the 
equilibria points as the roots of a polynomial of degree 4. These equations are 
solved and the critical colinear points are found to be unstable. 

In section 4 we use a continuity argument to prove the existence of an almost 
periodic orbit with an infinite number of collisions. The basic idea is to use the 
reciprocal of the distance of the center of mass of the two-satellite system to the 
center of force as a perturbation parameter e. For e = 0, the binary made by the 
satellites will be at an infinite distance of the center of force and the problem is 
integrable. Since equations are regularized, the Kepler equations for the binary are 
now represented by a ressonant harmonic oscillator (RHO). For small e, solutions 
of the perturbed problem are close to the solutions of the RHO. But solutions of 
the RHO that pass trough the origin are segments of straigth lines. Perturbing 
these solutions we expect that at least one of the perturbed solutions will preserve 
the feature of passing through the origin more then one time. This is the content 



of our main lemma (5.2). The symmetry of the Hamiltonian with respect to the 



origin guarantees the existence of the collision orbit. 

Acknowledgments. I am grateful to Jair Koiller for introducing me to the two- 
satellite problem. It's also a great pleasure to thank Hildeberto Cabral, Claudio 
Vidal and Marcelo Marchesin for many discussions and suggestions. I would like 
to thank CNPq-Brazil for financial support. 

2. The Model 

The Hamiltonian H of the planar problem of two bodies of mass mi and TO2 
attracted by a center of force of mass mg at the origin is 



2mi 2m2 \qi\ \q2\ m - q2\ 

where qi and 52 are the coordinates of the bodies of masses mi and m2 respectively, 
pi and p2 their conjugate momenta and G is the gravitational constant. We choose 
units such that G — I and set mo = 1. This Hamiltonian represents a Hill problem 
when TOi,TO2 << 1. We introduce a proportionality factor A G (0, cxd) such that 
1712 — Xrni. The Hamiltonian becomes 

- _ pI pI _ m]_ _ Ami _ \ml 

2mi 2Xmi \qi\ \q2\ 191-92!' 

Let w = dqi A dpi + dq2 A dp2 denote the standard symplectic 2-form. Let Xj^ be 
the Hamiltonian vector field generated by H. Consider the fiber scaling given by 

*(9i,92,Pi,P2) = {qi,q2,mipi,mip2). 



COLLISION ORBITS 



Xmi 
kl -92I 

w = rrii (dqi A dpi + dq2 A dp2) ■ 

Dividing Hamilton's equations ixgW = dH by mi we see that it suffices to study 
the Hamiltonian flow given by the Hamiltonian 

n? no 1 A Ami 

(1) H = i^ + ^ 



Under this scahng we have 






fr fp' ^ P' 


1 

kil 


A 

192 


and 







2 2A \qi\ \q2\ \qi-q2\ 

with standard symplectic 2-forni w = dqi A dpi + dq2 A dp2 ■ 
We introduce Jacobi variables p and r by 

91 = P - TTA ^' 

Here p represents the position of the center of mass of the two satellites and r 
represents their relative position vector. The Hamiltonian (|l|) becomes 



(2) H = 



Pp p; I A Ami 



2(1 + A) 2r \p-j^r\ \p 



A ^1 l^_L 1_ 



T+A' I If" ' 1 + X 



r 



where Pp and Pr are the momenta canonicaly conjugate to p and r respectively and 
r = j^ . Assuming -pj < ^-j^ we have the convergent expansions | Br | 

r-l 

1+A 



(3) T = -r-r y ^n(cos„, , , 



where P„ (x) is the n-th Legendre polynomial and 9 is the positively oriented angle 
between per. Hamiltonian ^) can be written as 

where 

A„-r"(l + (-!)" A), 

and A = 1 + A. 

The first parenthesis term of (Q) represents the Kepler problem described by the 
center of mass around the center of force and the second parenthesis term of fl) 
represents the Kepler problem of the two satellites around their center of mass. At 
this point we make the principal assumption of this work, namely, we assume that 
P = iPxTPy)j the vector representing the position of the center of mass of the two 
satellites, describes a circular keplerian orbit of radius \po\ around the center of 
force, i.e. p is a circular solution of p = ^rj^P yielding 

p= IpoI (cos(wt),sin(wt)), 

where uj — |/9o|^^- This is equivalent to consider the Circular Hill's Problem. By 
the second law of Kepler the energy of the center of mass is given by Ecm = — r^ , 
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and Hamiltonian (^) becomes 

(5) H^-^ + (§-^)'±P.icose)(m\,^. 

\po\ V2r \r\ J \po\ ^^ VIpoI/ 

We remark that this Hamihonian is time dependent since the angle 9 depends 
expUcitely on time. For future reference we write 

r.j, p^ + Ty py r.j, cos{ojt) + Vy sin(a;i) 

(b) cos(6') = — j-j— = —j . 

\r\\p\ \r\ 

Since energy of system (||) is not preserved we extend phase space from M^ to M^ 
by including the canonically conjugated pair {E,t). Our new Hamiltonian system 
is given by 

(7) |^ = -^-W+(#-Tr)-wS:r.r^.(cos^)(a)"A„, 
I w — du A dpu + dv /\ dpv + dE A di, 



where we must restrict our attention to the level set 7i = 0. Denoting the new time 

dt 
df 



by / it follows from Hamilton's equation ^ = 1. By choice we identify / and t. 



3. Regularization 
We regularize the collision between the two satellites. Writing p — {p.j.,,py) and 



r = {rx,ry) we write the Levi-Civita transformation |SS| 

,2 



(8) 



I X 



V 



Ty — 2UV, 



. Py^2wZ. 

Identifying r and p with the complex vectors r^+iry and Px+i Py respectively, we 
have that transformation (J^) can be written as r — {u + iv) and p — {w -\' i z) . 
This transformation halves angles and therefore takes the angle 9 between p and 
r to its half. Writing ^ = {u^v) and 7 = (tf,z) the transformation takes \r\ 
and \p\ to l^p and |7p respectively. Observe that under this transformation the 
curve p{t) — \pQ\{cos{wt),sa\{wt)) becomes "f{t) — \pq\'^ {cos{wt/2),sm.{wt/2)). 
Considering the lift of (||) to the cotangent bundle, Hamiltonian (|j) becomes 

T 1 „2 00 X 1^1 s 2n 

(9) 7i = -i?-^ + — (i^-Ami)-^$]P„(cos(0/2))^ A„. 

l7or W sr \^o\ t^i \m\J 

Observe that w = |po|~^ = l7o|~^ Regularization is achieved doing the time 
reparametrization given by 

where s denotes the new independent variable. This reparametrization can be 
performed considering the Hamiltonian 

(10) -H^\en 

in extended phase space with symplectic 2-form given by 
(11) w = du A dpu + dv A dpv + dt A dE. 
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Since the hypersurfaces {H = 0} and {7Y — 0} are equal it follows that the Hamil- 
tonian flow of (m at the level set 7i = is a reparametrization of the Hamiltonian 
flow of doh at the level set H = 0. M) and (M) yields 



We are interested on the flow of H at the level 0. We can eliminate the constant 
—Ami of the Hamiltonian by considering the level Xnii instead. The square of 
the reciprocal of the radius of the center of mass will be treated as a perturbation 
parameter. Writing e = -i — p- and doing the symplectic scaling p^ — > 2p^ , ^ — > ^/2 
we have 

(12) Hxm, ^ ^ ^ j{~Xe + E)\e ~ e' J2e-' P,,{cos{e/2)) r-SL\ A„, 

n=l ^ ^ 

where the subscript Xnii is a reminder that we must consider the level set Ti = Xmi. 

Lemma 3.1. Hamiltonian i\l^ ) is invariant with respect to the symmetry 5 : R^ — > 
M6 given by Si^,p^,E,t) = i-^,p^,E,t), i.e. n{^,p^,E,t) ='H{-^,p^,E,t) . 

Proof. Since |^| is invariant under 5 it sufflces to show that cos {6/2) is also invari- 
ant. But from (|) and (|) 

, ,^ iu^ -v'^)coB{wtl2) + 2uvsa\{wt/2) 
cos(0/2) = ^^ 

that is clearly invariant under S. D 

The equations of motion of ( |l2| ) can be written as 



(13) 



where 



^=-A|±S^ + g2yy. 



E^ „g2 ay. 



Mi! 

4 



F = ^e""ip„(cos(0/2))^ A„. 

n=l ^ '' 

From the second equation of ( [13| ) we have that E{s) — E{0) +0{e2)^ and we write 
for future reference that 

(14) j = _|25_4,,e>,a). 

4. Euler's Critical Points 

Euler's critical points are relative equilibria of the Hamiltonian system represent- 
ing colinear configurations. In what follows we compute the Euler's critical points 
for our model: In a colinear configuration the two masses will be moving forming 
a straigth line with the center of force. Therefore we must have 9 — or 6 = tt 
and 9 = 0. The first case represents the case where mi is between the center of 
force and m2 and the second case represents the case where m2 is in between. Since 
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m2 — \m\ we need to work only with the case = 0. We introduce polar coordi- 
nates (^ (/>) in the plane (u,w), where I is the radius and </> is the angle. Observe 
that (j) = TT— ""^"^ . Let pi and p^ be the conjugate momenta to I and (j) respectively. 
Hamiltonian ( [l2|) becomes 

(15) 



Ha™, = ||; + ^ - ^ (Ae + ^) ^' - e' E ^""'^" (cos(0/2)) r-] A„, 

n— 1 ^ ''^ 



Hamilton's equations are 
(16) 

( j — EL 
I pi 



Pi = i(Ae + i?)? + j?f + 2e2z^-^e"-ip„(cos(0/2)) (n + l) (4) A„, 



ri2' 



V) i?,P„(cos(0/2))sin(0/2)A„, 

(2 \ n+1 
V) 7^,P„(cos(0/2))sin(0/2)A„, 



, '- 4 1 

where 9/2 — tt — (j) — ^. We look for Euler's critical points of (iq). The condition 
^ = implies that (j) — wi/2. Using the second and last equations of (ll^) it follows 
that p^ = ™g . This is the value the angular momentum p^ must have in order 
to keep the colinear shape of the configuration. Since we are looking for colinear 
critical points in phase space (and not in extended phase space!) we must find a 
point for which I = 0, pi = 0, ptf, — and E = 0. For = i.e., for </> = 'wt/2 it 
follows from dlq) that the last two equalities are satisfied. We write E[) = £^(0). 
The first equality will be satisfied setting pi = 0. It remains to find I for which 
pi = 0. For 6* = we have that P„ (cos(|)) ~ -Pn(l) = 1- Recalling that w ^ e^ it 
follows that 

. dV 

(17) m = -^, 

where the potential function V is written as 
(18) 

Suming the series we obtain 

1 ^^ r. ^ ,2 e*r'^ e2;4r r4A + e;2r(l-A) 
(19) V{l)^~-{\e + Ea) P ' ' 




4 ^ ^ 16 4 1 16-e2Z4r2 

dV 
dl 



We want to find non-zero solutions of ^ = 0. Making u — P we only need to find 
non-zero solutions of ^ = 0. After some straightforwad computations we see that 
the zeros of this equation is given by the zeros of a polynomial of degree 4 in m. 
Therefore we have explicit analytic solutions for the Inner Euler's Critical points. 
The explicit solutions can be computed in an algebraic manipulator software but 
are very messy to be of any use. The critical points can be numerically computed 
using a Newton algorithm. We computed the values for diferent values of Eq. For 
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all of the critical points computed we have that ^;^ < indicating that the critical 
points are unstable. 

5. Collision Orbits 



Hamiltonian (12) is symmetric with respect to reflections about the origin. This 
simmetry will be used to prove existence of an almost periodic orbit that passes 
through the origin twice. Since the origin in the regularized plane represents a 
collision in the physical plane, this orbit represents a periodic orbit with an infinite 
number of collisions. The proof relies on the fact that for e = and Eq < 0, 
( [l2| ) represents a ressonant harmonic oscillator. Therefore, the projections on the 
{u, v) plane of trajectories leaving the origin (ejection trajectories) are segments of 
straight lines. We expect that when e is small, some trajectories of the perturbed 
system will preserve the feature of passing through the origin at least two times. 
The symmetry with respect to the origin will then imply that this orbit will cross 
the origin an infinite number of times. 

We set TOi = /x and assume that rrii, Ami << 1 and Eq < 0. We write ^ — 0{e). 
More explicitly we will assume that there is a constant p, such that /x(e) = efi. We 
will prove the following theorem 

Theorem 5.1. If e is small enough, there is an almost periodic orbit of ( fZ^ j that 
passes through the origin. 



Proof. For easy reference we write Hamiltonian (|12|) 

(20) H,^^pl-^CXe + E)\^\'-e^V, 

where k — Xfi. We restrict our study to solutions that at time s = leave the 
origin. 

Definition We call ejection solutions, solutions of the Hamiltonian system with 
Hamiltonian (pfl) and symplectic 2-form (O) that at time s = leave the origin, 
i.e., solutions with ^(0) = (0,0). 

Remark: The set of ejection trajectories is parametrized by a circle. In fact, 
for these trajectories ( |2(]| ) implies that |p^(0)p — en and we can write Pu(0) = 
y/eKCOs{a) and PviO) — \/f^sin{a) for a G [0, 27r) . 

For a small e and £'o < solutions of ( pO| ) will be close to the solutions of a 
ressonant harmonic oscillator with period and amplitude given by 



27rV2r J , „ I en 

In = — • and A 



VW\ V l^ol 

respectively. Let r/ denote an initial condition for an ejection trajectory, i.e. rj ~ 
(0, 0,Pu,Pv, E{0), t(0)). Fixing £'(0) and t{0) we have that 77 is uniquely determined 
by a. We write the set of ejection trajectories at time s and parameter e as 

0'(a) = {<ia),v'^{a),Putia),pytia),E'^{a),tl{a)). 

Then <f>°{a) = ?? == (0, 0, y/n^coa{a), y/el^sin{a),E{0),t{0)). 

The following lemma shows that, for e small enough, an ejection trajectory with 
angle a pointing to the right half plane will cross the v axis transversaly at some 
time T = T{a). 
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Lemma 5.2. Let Eg < 0. Let 6 be a positive real number with 6 << 1. Let a £ I, 

where I — [— -^ + 5, ^ — 5] . Lf e is small enough, there exist times t = T(e, a) such 
that 

(21) ^To<T<^To, 

and uj(a) = . Moreover v^{a) is a continuous function of a. 

Proof. Wc first prove existence of t. Let m ~ min{| cos(a)|} = | cos(— — (5)|- We 
write 

(22) (bt=ros, + er 

wfiere (t)osc i^ tlie liarmonic oscillator flow given by 

4'lsc(.^) — (A cos(a)sin(aJs), A sin(a)sin(ijjs),— Aijjcos(a) cos(a;s), 
-Aujsin{a) cos{ujs),E{0),t{0)), 

where w = •\/2r • F^'^i^ (111) "01 satisfies 

En 0? 

Let t == i To and f = I Tq. Then 

(23) Uj{a) = y^Tfi cos{a) + e'il!-{a); 

where "0^ is the u component of "0*. Let J — [0,e] where e is a positive small 
number. Let M = max I max |i/)-(a)| I . We can find e small enough such that 

(24) e^ k"^ m> eM, 
this implies that 

|v^cos(a)| >e|V'|(a)|. 

Therefore for e satisfying (gj), (g^) implies that ul has the same sign as ^/Tk cos(q;), 
i.e., u* is positive. Analogously we show that u* is negative. Therefore by continuity 
of the flow, there exists times t^, for which uj — and satisfying (21). To prove 
continuity it suffices to prove that at time s — t the projection of the flow on the 
{u, v) plane intersects the v axis transversaly, i.e., it suffices to show that ul{a) ^ 0. 
But 

ul{a) — — ^wsin(Q;) cos(u;r) + eip^. 

For e small enough it follows from ( p4|) , using the same argument as in the existence 
part, that u1{a) ^0. D 

For notational simplicity we will write in the next lemma u'l — u{s) and w| = ^(s). 

Lemma 5.3. Let Eq < 0. Then for e small enough there exists a G / such that 
v{Tuia)) = 0. 



Proof. Let r as in lemma ( ]5.2D . By Taylor's Theorem it follows that there exists c, 
< c < T such that £_{t) = ^(0) + C(O) t+^t"^ + ^^r^, that we write as 

u{t) = u{0) + u{Q) r + ^T^+'^T^ 
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v{t) = v{0) + m r+^r'+'^T^ 

For an ejection trajectory we have that u{0) = v{0) — and equations ( |l^ ) imply 
that u(0) = v{0) = 0. Therefore 

u(t) = \JtK cos(a) T H r , 

6 

w(r) = V^'^smfajr H — t , 

6 

with < c < r. r satisfies ( |2l| ) and by definition M(r(a;)) = 0. Solving the first 
equation for t and substituting on the second equation we obtain 



(25) w(T„(a)) = 7777^ (sin(a)M (c) - cos(a)'t''(c)) . 

m(c) 

Now we estimate | ii (c) — v (c) | . From the Taylor expansions we have that 

6 6 

l'w(c) - 'v{c)\ < -^ \^/eK{cos{a) - sin(a))| + -^ \u{t) - v{t)\. 



Using (^ we obtain 



6V^,^ , , 6|go|J 



\uic) - vic)\ < -^ \Eo\ + ^^ Hr)\. 



But, from (E2|) we have that v{t) ~ t^-t sin(a) cos(ws) + 0{e ) giving that 

\uic)-v{c)\<-^V^\Eo\ + ^^^ + Oie'). 
T^^ 7r2v2 \Eq\2 

Therefore there exists a constant K such that 

\u{c) - v{c)\ <K\Eo\ ^ + 0(e3). 

We can then write 

■w(c) = w(c)+0(e5), 

and ( psi) becomes 



v{Tu{a)) = ^1^ (sin(a) - cos(a) + ©(e^ ) j . 

Thus for e small enough we can find ai, a2 G / such that w(t(q;i)) > and 
w(r(a2)) < 0. By the continuity part of lemma ( ^.2| ) there is a € /, such that 
v{T(a)) = 0. D 

Therefore we proved that for e small enough there is an angle a such that the 
ejection trajectory with this angle will pass through the origin in some future time 
T(a). But by lemma (3.1) the system is symmetric with reflections about the origin. 
Therefore, the orbit will pass an infinite number of times around the origin proving 
the theorem. D 
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